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Abstract 


Ocean  wave  propagation  causes  random  change  of  ocean  surface  slope  and  in 
turn  affects  the  underwater  bomb  trajectory  deviation  (r)  through  water  column.  This 
trajectory  deviation  is  crucial  for  the  clearance  of  obstacles  such  as  sea  mines  or  maritime 
improvised  explosive  device  (IED)  in  coastal  oceans  using  bombs.  A  nonlinear  six 
degrees  of  freedom  (6-DOF)  model  has  been  recently  developed  and  verified  at  the  Naval 
Postgraduate  School  with  various  surface  impact  speeds  and  surface  slopes  as  model 
inputs.  The  surface  slope  (5)  randomly  changes  between  0  and  nil  with  a  probability 
density  function  (PDF)  p{ s),  or  called  the  5-PDF.  After  5  is  discretized  into  /  intervals  by 
5 1,  52,  ...,  Si ,  ...  5/+1,  the  6-DOF  model  is  integrated  with  a  given  surface  impact  speed 
(v0)  and  each  slope  s,  to  get  bomb  trajectory  deviation  r  at  depth  (h)  as  a  model  output. 

The  calculated  series  of  { r  }  is  re-arranged  into  monotonically  increasing  order  (r,} .  The 

bomb  trajectory  deviation  r  within  (rh  rJ+\)  may  correspond  to  one  interval  or  several 
intervals  of  5.  The  probability  of  r  falling  into  (17,  rr  1 )  can  be  obtained  from  the 
probability  of  5,  and  in  turn  the  PDF  of  r,  or  called  the  r-PDF.  Change  of  the  r-PDF 
versus  features  of  the  5-PDF,  water  depth,  and  surface  impact  speed  is  also  investigated. 

Keywords:  3D  underwater  bomb  trajectory  model,  probability  density  function,  bomb 
trajectory  deviation,  stochastic  ocean  surface  slope,  STRIKE35 
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1.  Introduction 


Movement  of  a  fast-moving  rigid  body  such  as  a  bomb  through  water  column  has 
been  studied  recently  [1-3].  These  studies  have  been  motivated  by  a  new  concept  of  using 
the  Joint  Direct  Attack  Munition  (JDAM,  i.e.,  ‘smart’  bomb  guided  to  its  target  by  an 
integrated  inertial  guidance  system  coupled  with  a  global  positioning  system)  Assault 
Breaching  System  (JABS)  for  mine/maritime  improvised  explosive  device  (IED) 
clearance,  in  order  to  reduce  the  risk  to  personnel  and  to  decrease  the  sweep  timeline 
without  sacrificing  effectiveness  (Fig.  1).  Underwater  bomb  trajectory  depends  largely 
on  the  surface  impact  speed  and  angle.  When  the  surface  impact  of  high-speed  rigid  body 
such  as  scaled  MK-84  warhead  is  normal  or  near  normal  to  the  flat  water  surface,  four 
types  of  trajectories  have  been  identified  from  experimental  and  numerical  modeling 
results  [4]  depending  the  characteristics  of  the  warheads:  with  tail  section  and  four  fins 
(Type-1),  with  tail  section  and  two  fins  (Type-11),  with  tail  section  and  no  fin  (Type-111), 
and  with  no  tail  section  (Type-IV).  Type-1  trajectories  are  quite  stable  downward 
without  oscillation  and  tumbling  no  matter  the  water  entry  velocity  is  high  or  low.  Type-2 
and  Type-3  trajectories  are  first  downward,  then  makingl80°  turn  (upward),  and  travel 
toward  the  surface.  Type-IV  trajectories  are  at  first  downward  with  little  horizontal  drift 
and  then  tumbling  downward  with  large  horizontal  drift. 

The  horizontal  distance  (r)  (or  called  trajectory  deviation)  between  surface  impact 
point  and  the  bomb  location  varies  with  depth  in  different  types  of  trajectories.  This 
parameter  draws  attention  to  the  naval  research  due  to  the  threat  of  mine  and  maritime 
IED.  Prediction  of  trajectory  deviation  of  an  underwater  bomb  contributes  to  the  bomb 
breaching  for  mine  and  maritime  IED  clearance  in  surf  and  very  shallow  water  zones 
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with  depth  shallower  than  12.2  m  (i.e.,  40  ft),  shallow  water  zones  (12.2  -  91.4  m,  i.e., 
40-300  ft),  and  deep  zones  (deeper  than  91.4  m,  i.e.,  300  ft)  according  to  U.S.  Navy’s 
standards.  The  bombs’  trajectory  drift  is  required  to  satisfy  the  condition,  r  <  2.1  m,  for 
the  validity  of  mine  clearance  using  bombs  [5], 

In  coastal  oceans,  waves  form  when  the  water  surface  is  disturbed,  for  example, 
by  wind  or  gravitational  forces.  During  such  disturbances  energy  and  momentum  are 
transferred  to  the  water  mass  and  sea-state  is  changed.  For  very  shallow  and  shallow 
water  regions,  the  bottom  topography  affects  the  waves  dramatically  and  causes  a 
significant  change  in  surface  slope.  When  bomb  strikes  on  the  wavy  ocean  surface,  a 
scientific  problem  arises:  How  does  randomly  changing  ocean-surface  slope  affect  the 
underwater  bomb  trajectory  and  orientation?  Or  what  is  the  probability  density  function 
of  the  underwater  bomb  trajectory  deviation  due  to  random  sea  surface  slope?  The  major 
task  of  this  paper  is  to  answer  these  questions.  The  effect  of  surface  slope  on  the 
underwater  bomb  trajectory  is  presented  in  Section  2.  Stochastic  features  of  the  sea  slope 
are  simply  described  in  Section  3.  A  recently  developed  six  degrees  of  freedom  (6-DOF) 
model  at  the  Naval  Postgraduate  School  for  predicting  underwater  bomb  location  and 
trajectory  is  depicted  in  Section  4.  Ensemble  6-DOF  modeling  to  get  PDF  of  trajectory 
deviation  d  from  the  stochastically  changing  sea  surface  slope  is  described  in  Section  5. 
The  conclusions  are  presented  in  Section  6. 

2.  Effect  of  Ocean  Surface  Slope  on  Underwater  Bomb  Trajectory 
Let  /j  be  the  inclination  angle  of  the  ocean  surface;  and  (f)  be  the  bomb  impact 
angle  relative  to  the  normal  direction  of  the  ocean  surface  (Fig.  2).  For  a  flat  surface 
(no  waves), 
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(la) 


//  =  0 

For  90°  bomb  striking  (vertically  downward), 

0  =  0  (lb) 

With  ocean  wave  propagation,  //  can  be  treated  as  an  averaged  value  in  a  wave  period; 
and  corresponding  averaged  slope  in  a  wave  period  (5)  is  given  by 

5*  =  tan  /u .  (2) 

The  ocean  waves  may  cause  evident  slant  of  the  ocean  surface  with  /u  «  55°  (Kinsman, 
1965),  which  affects  the  underwater  bomb  trajectory,  orientation,  and  horizontal  drift  (r) 
(Fig.  3).  The  differential  effects  depend  on  which  part  of  the  wave  is  impacted  by  the 
bomb  (i.e.,  different  sea  slopes).  Obviously,  such  a  wave  effect  can  be  investigated  by  a 
6-DOF  model  with  a  sloping  surface  (i.e.,  //  changing  with  time)  and  non-nonnal 
impact  angle  (i.e.,  0  0  0). 

Besides,  the  surface  slope  also  affects  the  tail  separation  due  to  the  bomb  and 
cavity  orientations  and  the  air-cavity  geometry.  This  is  because  the  air  cavitation  or 
supercavitation  is  usually  generated  after  the  bomb  enters  the  water  surface  [7].  The 
cavity  is  usually  oriented  in  the  same  direction  of  the  bomb  velocity  with  its  geometry 
simply  represented  by  a  cone  with  the  angle  (7).  The  bomb  orientation  relative  to  the 
cavity  is  represented  by  the  angle  between  the  bomb  main  axis  and  velocity  (/?).  The 
condition  for  bomb  not  hitting  the  cavity  wall  is  given  by  (Fig.  4a) 

P  <  Y  •  (3a) 
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Violation  of  the  condition  (3)  may  cause  the  tail  separation  (bomb  hitting  the  cavity  wall), 
as  shown  in  Fig.  4b.  Ocean  waves  not  only  affect  the  bomb  trajectory  and  orientation  but 
also  change  the  cavity  orientation,  which  may  cause 

P>Y ,  (3b) 

i.e.,  the  bomb  may  hit  the  cavity  wall  and  cause  the  tail  separation  (Fig.  4). 

3.  PDF  of  Ocean  Surface  Slope 

Wave  height  and  wave  period  are  approximately  independent  of  each  other  for 
either  wind  waves  or  swells,  but  not  for  mixed  waves.  From  mixed  wave  records,  Gooda 
[8]  found  that  there  is  a  strong  correlation  between  wave  height  and  wave  period.  In  fact, 
the  correlation  is  mainly  caused  by  the  two  or  more  groups  of  notable  waves  with 
different  characteristic  wave  heights  and  periods  in  the  mixed  waves.  With  the 
independent  assumption  between  wave  amplitude  and  wave  period  (or  wavelength),  the 
PDF  of  averaged  wave  slope  5  scaled  by  its  standard  deviation  a  (the  real  slope  is  s*  =  so) 
is  obtained  from  the  PDF  of  wave  length  and  PDF  of  wave  amplitude  [9], 

-(«+ 2)12 

(4) 

where  n  is  the  peakedness  coefficient  which  is  determined  by  both  the  spectral  width  of 

the  gravity  waves,  and  the  ratio  between  the  gravity  wave  mean-square  slope  and  the 

detectable  short  wave  mean-square  slope.  Generally  speaking,  the  peakedness  of  slopes 

is  generated  by  nonlinear  wave-wave  interactions  in  the  range  of  gravity  waves;  and  the 

skewness  of  slopes  is  generated  by  nonlinear  coupling  between  the  short  waves  and  the 

underlying  long  waves.  For  n  =  2,  the  PDF  of  the  wavelength  corresponds  to  the  Rayleigh 

distribution.  For  n  =  10,  the  PDF  in  (4)  fits  the  Gram  Charlier  distribution  [10],  very 
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well  in  the  range  of  small  slopes.  As  «— > oo,  the  PDF  of  the  wavelength  tends  to  the 
Gaussian  distribution  [9].  Fig.  6  shows  four  typical  surface-slope  characteristics:  (a)  n  = 
2,  (b)  n  =  4,  (c)  n  =  10,  and  n  =  100.  It  is  seen  that  There  is  almost  no  difference  in  PDF 
between  n  =  10  and  n  =  100. 

4.  A  6-DOF  Model  (STRIKE35) 

Recently,  a  6-DOF  model  has  been  developed  at  the  Naval  Postgraduate  School 
for  predicting  underwater  bomb  location  and  trajectory.  It  contains  three  parts: 
momentum  equation,  moment  of  momentum  equation,  and  semi-empirical  formulas  for 
drag,  lift,  and  torque  coefficients  [11-13].  The  momentum  equation  of  a  rigid  body  is 
given  by 

'"*  =  F«+F„  +  F,+r,,  (5) 

where  m  is  the  mass  of  the  rigid  body,  u  is  the  translation  velocity  of  the  center  of  mass, 

¥g  =  -mgk,  ¥h  =  pTlgk,  (6) 

are  the  gravity  and  buoyancy  force;  II  is  the  volume  of  the  rigid  body;  k  is  the  unit 
vector  in  the  vertical  direction  (positive  upward):  and  g  is  the  gravitational  acceleration. 
Frf  is  the  drag  force;  and  F/  is  the  lift  force. 

The  moment  of  momentum  equation  is  given  by 

J#^  =  -^ex(pngk)  +  M/,,  (7) 

where  Q  is  the  rigid-body’s  angular  velocity  vector;  cr  is  the  distance  between  center 
of  volume  (ov)  and  center  of  mass  (om),  which  has  a  positive  (negative)  value  when  the 


7 


Launch  Internet  Explorer  Browser. Ink 


direction  from  ov  to  om  is  the  same  (opposite)  as  the  unit  vector  e;  Mh  is  the 
hydrodynamic  torque  due  to  the  drag/lift  forces;  and  J  is  the  gyration  tensor. 

The  drag/lift/torque  coefficients  should  be  given  before  running  the  6-DOF  model. 
These  coefficients  depend  on  various  physical  processes  such  as  water  surface 
penetration,  super-cavitation,  and  bubble  dynamics.  A  diagnostic-photographic  method 
has  been  developed  [4]  to  get  semi-empirical  formulae  for  calculating  the  drag/lift/torque 
coefficients  for  underwater  bombs  with  dependence  on  the  Reynolds  number  (Re),  angle 
of  attack  (a),  and  rotation  rate  along  the  bomb’s  major  axis  (Q)  [4], 


Here,  Re*  =  1.8  x  1 07,  is  the  critical  Reynolds  number,  and 

i 

d  =  [n:2'2  -(^-|;r-2<2|)~~j2'2sign(;r-2a),  (11) 

(2  aV'8  (2a  Y'7 

6x=k  —  ,  d2  =  In  — -1  .  (12) 

V  ft  J  V  n  v 

5.  PDF  of  Bomb’s  Horizontal  Drift 
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Let  the  bomb  be  dropped  in  the  vertical  direction  to  the  slanted  sea  surface 
characterized  by  an  averaged  slope  (5*  =  a.s)  in  a  wave  period,  here  5*  =  tan  p  (see  Fig. 
1).  Consider  a  5-time  of  s*  value  as  the  interval  [0,  5s*]  for  the  change  of  the  surface 
slope.  This  interval  [0,  5s*]  is  divided  into  /  equal  sub-intervals, 


5/v  * 

CJS,=—,  /  =  0,1,2,...,/. 


(13) 


with  the  corresponding  inclination, 


//  =  arc  tan  (cry)  =  arctan 


5  is '' 


i  =  0,1,2,...,/. 


(14) 


For  a  given  parameter  n  in  the  5-PDF,  the  probability  for  5*  taking  values  between  asy  1 
and  osj  is  calculated  by 

P  =  Prob(y  <  s  <  y+1)  =  f  +'  p(s)ds  .  (15) 

J  Sj 

The  6-DOF  model  is  integrated  /  times  (called  ensemble  integration)  from  the 
surface  impact  speed  ( V)  and  various  p,  values  to  get  the  bomb  horizontal  drift  r  (i  =  0, 

1,  ...,/)  at  depth  z  =  -H.  The  series  {  rt ,  i  =  0,  1,  . . .,  /}  might  not  be  in  monotonically 

increasing  or  decreasing  order.  Therefore,  it  is  reorganized  into  monotonically  increasing 
order  {y  J  =  0,  1,  ...,  J }  with  J  <  I.  The  inequality  is  due  to  an  interval  [rn  rj+{\  of  the 
horizontal  drift  corresponding  to  m  intervals  {[.v,i,  sn+1],  [s/2,  5,2+ 1],  ...,  [sim,  sim+ 1  ] }  of 
the  surface  slope  (Fig.  7)  .  The  probability  for  the  bomb’s  horizontal  drift  r  taking  values 
between  r7-  and  rJ+  \  is  calculated  by 

Q  =  Prob(r,  <r<  r.+] )  =  f  1+1  p(s)ds  +  f  '  p(s)ds  + ...  +  [  p(s)ds  .  (16) 

^  ^  "  si\  J  S' 2  J  Sjn 
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The  probability  density  between  r7-  and  rJ+\  is  calculated  by 


Pj 


(17) 


From  pj,  we  can  obtain  the  PDF  of  r,  or  called  the  r-PDF. 

6.  Sensitivity  Studies 

6.1  Dependence  of  the  r-PDF  on  Depth 

Dependence  of  r-PDF  on  depth  can  be  identified  from  the  ensemble  integration  (/ 
=  100)  of  the  6-DOF  model  with  given  bomb’s  surface  impact  speed  (V=  300  m/s),  5*  = 
0.2  (i.e.,  a  =  0.2),  and  n  =  2  (i.e.,  large  peakedness  in  the  5-PDF).  The  calculated  r-PDF 
(Fig.  8)  is  positively  skewed  for  shallow  depth  ( H  =  12.2  m,  i.e.,  40  ft),  reduces  the 
skewness  as  depth  increases  to  50  m,  becomes  negatively  skewed  as  the  depth  exceeding 
91.4  m  (i.e.,  300  ft).  The  negative  skewness  strengthens  as  depth  deeper  than  91.4  m. 
The  horizontal  axis  in  all  the  panels  Fig.  8  is  the  non-dimensional  horizontal  drift  r!H. 
The  median  (50  percentile  <70.5)  of  the  horizontal  drift  (r)  is  0.16  m  at  the  depth  z  =  -12.2 
m,  1.7  m  at  z  =  -50  m,  5.4  m  at  z  =  -91.4  m  (300  ft),  18.0  m  at  z  =  -150  m,  34.0  m  at  z  =  - 
200  m,  and  52.5  m  at  z  =  -250  m  (Table  1).  Here  z  is  the  vertical  coordinates  with  z  =  0 
corresponding  to  the  water  surface.  Thus,  down  to  the  depth  of  50  m,  the  median  value  of 
the  horizontal  drift  is  always  less  than  the  Navy’s  criterion,  i.e.,  2.1  m.  The  95  percentile 
(<70.95)  of  the  horizontal  drift  (r)  represents  a  reasonable  estimation  (with  95%  of 
confidence)  of  the  distance  between  bomb  and  mine/maritime  IED  when  the  bomb 
maneuvering  in  the  water  column.  If  this  value  is  smaller  than  2.1  m,  according  to  the 
Navy’s  standard,  the  bomb  will  effectively  ‘kill’  the  mine/maritime  IED.  It  is  0.32  m  at 
the  depth  z  =  -12.2  m,  2.8  m  at  z  =  -50  m,  7.86  m  at  z  =  -91.4  m  (300  ft),  22.5  m  at  z  = 


10 


Launch  Internet  Explorer  Browser. Ink 


-150  m,  40.0  m  at  z  =  -200  m,  and  60.0  m  at  z  =  -250  m  (Table  2).  The  5  percentile  (<70.05) 
of  the  horizontal  drift  (r)  represents  the  minimum  distance  (likely)  between  bomb  and 
mine/maritime  IED  when  the  bomb  maneuvering  in  the  water  column.  It  is  0.13  m  at  the 
depth  z  =  -12.2  m,  0.6  m  at  z  =  -50  m,  5.48  m  at  z  =  -91.4  m  (300  ft),  10.5  m  at  z  =  -150 
m,  24.0  m  at  z  =  -200  m,  and  40.0  m  at  z  =  -250  m  (Table  3). 

6.2.  Dependence  of  the  r-PDF  on  the  Peakedness  of  the  5-PDF 

Keeping  all  the  initial  conditions  in  running  the  6-DOF  model  the  same  as 
described  in  subsection  6.1  except  changing  the  parameter  n  of  the  s-PDF  from  2  to  100 
(small  peakedness),  the  ensemble  integration  of  the  6-DOF  model  shows  the  following 
results.  The  calculated  r-PDF  (Fig.  9)  is  almost  zero  skewness  for  shallow  depths  (//  = 
12.2  m,  50  m),  becomes  negatively  skewed  as  the  depth  of  91.4  m  (i.e.,  300  ft).  The 
negative  skewness  strengthens  as  depth  deeper  than  91.4  m.  Comparing  between  Fig.  9 
and  Fig.  8,  we  may  find  that  the  negative  skewness  of  r-PDF  increases  as  n  increases. 
The  median,  <70.95,  and  <70.05  of  the  horizontal  drift  (r)  do  not  change  too  much  as  n 
increases  from  2  to  100  (Tables  1-3). 

6.3.  Dependence  of  the  r-PDF  on  the  Averaged  Surface  slope  a 

Keeping  all  the  initial  conditions  in  running  the  6-DOF  model  the  same  as 
described  in  subsection  6.1  except  increasing  the  averaged  surface  slope  a  from  0.2  to  1, 
the  calculated  r-PDF  (Fig.  10)  is  negatively  skewed  at  all  depths;  and  the  negative 
skewness  enhances  as  the  depth  increases.  Comparing  between  Fig.  10  and  Fig.  8,  we 
may  find  that  the  negative  skewness  of  r-PDF  increases  as  a  increases.  The  median,  <70.95, 
and  <70.05  of  the  horizontal  drift  (r)  increase  drastically  as  a  increases  from  0.2  to  1.0 
(Table3  1-3).  For  example,  <70.95  is  0.54  m  at  depth  z  =  -12.2  m,  4.0  m  at  z  =  -50  m,  10.05 
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m  at  z  =  -91.4  m  (300  ft),  25.5  m  at  z  =  -  150  m,  46.0  m  at  z  =  -200  m,  and  67.5  m  at  z  = 
-250  m  (Table  2). 

6.4.  Dependence  of  the  r-PDF  on  the  Surface  Impact  Speed  V 

Keeping  all  the  initial  conditions  in  running  the  6-DOF  model  the  same  as 
described  in  subsection  6.1  except  decreasing  the  surface  impact  speed  V  from  300  m/s  to 
200  m/s,  the  calculated  r-PDF  (Fig.  1 1)  is  quite  comparable  to  case  with  the  impact  speed 
of  300  m/s  (Fig.  8)  such  as  positive  skewness  for  shallow  depth  (H  =  12.2  m,  i.e.,  40  ft), 
weaker  skewness  as  depth  increasing  to  50  m,  negative  skewness  as  the  depth  exceeding 
91.4  m  (i.e.,  300  ft).  Comparing  between  Figs.  11  and  8,  reduction  of  surface  impact 
speed  leads  to  the  increase  of  the  peakedness  of  the  r-PDF.  The  median,  q0_ 95,  and  <70.05  of 
the  horizontal  drift  (r)  are  usually  higher  for  V  =  200  m/s  than  that  for  V  =  300  m/s  except 
for  the  very  shallow  water  depth  (z  =  -12.2  m)  where  <70.95,  and  <70.05  are  lower  for  V  = 
200  m/s  (0.17  m,  0.04  m)  than  that  for  V=  300  m/s  (0.32  m,  0.13  m)  (Tables  1-3). 

7.  Conclusions 

The  PDF  of  the  horizontal  drift  of  underwater  bomb  trajectory  (i.e.,  r-PDF)  due  to 
stochastic  ocean  surface  slope  is  obtained  through  ensemble  integration  of  the  6-DOF 
model  recently  developed  at  the  Naval  Postgraduate  School.  For  a  bomb  dropping  in  the 
vertical  direction  to  a  slanted  sea  surface,  the  input  parameters  of  the  6-DOF  model  are 
the  bomb’s  surface  impact  speed  (V),  and  surface  slope.  The  surface  slope  is  a  random 
variable  depending  on  two  parameters:  (a)  averaged  slope  within  a  wave  period  (a),  and 
(b)  peakedness  of  the  5-PDF  (n).  The  s-PDF  is  discretized  into  /  intervals  (in  this  paper,  / 
=  100).  For  given  values  of  (V,  a,  n),  the  6-DOF  model  is  integrated  /  times  with 
different  values  of  the  surface  slope  from  the  5-PDF  to  obtained  /  values  of  the  horizontal 
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drift  at  various  depth.  The  r-PDF  is  then  constructed  from  these  r  values.  The  r-PDF  has 
the  following  features: 

(1)  The  r-PDF  varies  with  depth.  Usually,  the  r-PDF  is  positively  skewed  for  very 
shallow  water  (77=  12.2  m,  i.e.,  40  ft),  and  negatively  skewed  down  below.  Increase  of 
the  peakedness  parameter  of  the  5-PDF  (n)  or  the  averaged  surface  slope  in  a  wave  period 
(a)  reduces  the  positive  skewness  at  the  very  shallow  water  and  enhances  the  negative 
skewness.  Decrease  of  the  bomb’s  surface  impact  speed  (V)  enhances  the  peakedness  of 
the  r-PDF.  Three  measures  were  calculated  (<70.05,  <70.5,  and  <70.95)  from  the  r-PDF. 

(2)  The  values  of  <70.95  are  small  for  all  cases  at  a  very  shallow  depth  (z  =  -12.2  m, 

1. e.,  40  ft)  with  a  maximum  value  of  0.54  m  for  the  initial  conditions  of  (V=  300  m/s,  n  = 

2,  a  =  1.0).  This  value  (0.54  m)  is  much  smaller  than  the  critical  value  of  2.1  m  for 
effectively  ‘killing’  the  mine/maritime.  This  may  prove  that  the  Joint  Direct  Attack 
Munition  (JDAM)  Assault  Breaching  System  (JABS)  is  effective  to  clear  mines  and  light 
obstacles  in  very  shallow  water  (depth  up  to  12.2  m,  i.e.,  40  ft). 

(3)  The  values  of  <70.95  are  all  larger  than  2.1  m  when  the  depth  deeper  than  50  m. 
This  indicates  that  to  extend  the  JABS  from  very  shallow  water  (12.2  m  depth)  to  shallow 
water  (12.2  m  -  91.4  m)  needs  more  studies. 
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Table  1.  The  median  horizontal  drift  (unit:  m)  of  an  underwater  bomb  at  various  depths 
obtained  from  ensemble  integration  of  the  6-DOF  model  with  various  input  parameters. 


Depth  (m) 

Case  1: 

V=  300  m/s 
n  =  2 
<7  =  0.2 

Case  2: 

V  =  300  m/s 

ft  =  100 

<7  =  0.2 

Case  3: 

V=  300 
m/s 

ft  =  2 

<7=  1.0 

Case  4: 

V  =  200  m/s 

ft  =  2 

<7=0.2 

12.2 

0.16 

0.16 

0.37 

0.17 

50.0 

1.7 

1.8 

3.1 

2.5 

91.4 

5.4 

5.7 

8.6 

8.9 

150.0 

18.0 

18.0 

22.5 

25.5 

200.0 

34.0 

34.0 

42.0 

44.0 

250.0 

52.5 

55.0 

62.5 

65.0 

Table  2.  The  values  of  <70.95  for  the  horizontal  drift  (unit:  m)  of  an  underwater  bomb  at 
various  depths  obtained  from  ensemble  integration  of  the  6-DOF  model  with  various 
input  parameters. 


Depth  (m) 

Case  1: 

V  =  300  m/s 

ft  =  2 

<7=0.2 

Case  2: 

V=  300 

m/s 

ft  =  100 

<7=0.2 

Case  3: 

V=  300 

m/s 

ft  =  2 

<7=  1.0 

Case  4: 

V=  200 
m/s 

ft  =  2 

<7  =  0.2 

12.2 

0.32 

0.27 

0.54 

0.17 

50.0 

2.8 

2.55 

4.0 

3.6 

91.4 

7.86 

7.40 

10.05 

10.97 

150.0 

22.5 

21.0 

25.5 

28.5 

200.0 

40.0 

38.0 

46.0 

48.0 

250.0 

60.0 

60.0 

67.5 

70.0 
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Table  3.  The  values  of  <70.05  for  the  horizontal  drift  (unit:  m)  of  an  underwater  bomb  at 
various  depths  obtained  from  ensemble  integration  of  the  6-DOF  model  with  various 
input  parameters. 


Depth  (m) 

Case  1: 

V  =  300  m/s 
n  =  2 
<7  =  0.2 

Case  2: 

V=  300 
m/s 

n  =  100 
<7  =  0.2 

Case  3: 

V=  300 
m/s 
n  =  2 
<7=  1.0 

Case  4: 

V=  200 
m/s 
n  =  2 
<7=0.2 

12.2 

0.13 

0.05 

0.15 

0.04 

50.0 

0.6 

0.80 

1.8 

1.05 

91.4 

5.48 

7.40 

5.76 

5.30 

150.0 

10.5 

12.45 

18.0 

18.0 

200.0 

24.0 

26.0 

34.0 

32.0 

250.0 

40.0 

45.0 

55.0 

55.0 
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BOMB  FALL  LINE 


Fig.  1 .  The  concept  of  airborne  sea  mine/maritime  IED  clearance. 
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Fig.  2.  Ocean  surface  inclination  angle  ( // )  and  bomb  impact  angle  ((f) )  relative  to  the 
nonnal  direction  of  the  surface. 


19 


Launch  Internet  Explorer  Browser. Ink 


I 


Fig.  3.  Dependence  of  underwater  bomb  trajectory,  orientation,  and  horizontal  deviation 
(, r )  on  the  ocean  surface  slope  or  on  different  locations  of  the  waves. 
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(b) 


(a) 


\  \ 

M  M 

Fig.  4.  Air  cavity  with  (a)  /3  <  y  (tail  section  not  hitting  the  cavity  wall,  and  (b)  with 
P  =  y  (tail  section  hitting  the  cavity  wall). 


21 


Launch  Internet  Explorer  Browser. Ink 


Fig.  5.  Wave  effect  on  the  air  cavity  orientation  which  may  cause  >  y  (Tail  section 
hitting  the  cavity  wall). 
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Probability  Density 


r/H 


r/H 


Fig.  11.  Probability  distribution  of  the  bomb’s  horizontal  drift  (scaled  by  the  depth)  r/H 
with  n  =  2,  a  =  0.2,  V  =  200  m/s  for  various  depth:  (a)  12.2  m  (i.e.  40  ft),  (b)  50  m,  (c) 
91.4  m  (i.e.,  300  ft),  (d)  150  m,  (e)  200  m,  and  (f)  250  m. 
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